F. Shirzaditabar; B. Oskooi; M. Bastani
Abstract
Airborne electromagnetic data can be obtained in time and frequency domains. In time domain, data consist of the time derivative of the secondary field while in frequency domain it is the absolute value of the secondary field. Here, we have used Occam’s inversion method on airborne electromagnetic ...
Read More
Airborne electromagnetic data can be obtained in time and frequency domains. In time domain, data consist of the time derivative of the secondary field while in frequency domain it is the absolute value of the secondary field. Here, we have used Occam’s inversion method on airborne electromagnetic data to resolve the resistivity of each layer in an assumed 1D layered earth. In the current modeling the number of model parameters is much more than the number of data. So, by choosing an appropriate objective function to be minimized, we have used the least squares method such that model parameters reached to a minimum structure and data are sufficiently fitted. One of the advantages of this method is that it doesn’t need a primary model and converges to the final model for every starting model. Applying this method on synthetic data provides reasonable results. The method is applied on real airborne electromagnetic data from Kalat-e-Reshm area in the east of Semnan province in Iran. The recovered maps and sections prove the existence of an arc like structure which has outcrop at the surface and is shown in the geological map. Besides, it reveals that the structure is considerably thick due to the alluvium deposits with 10 m thickness covering the area with a dip towards south.
B. Oskooi; G.A. Fanaee-Kheirabad
Abstract
Geothermal energy is playing a large role as an alternative energy source for both electricity generation and space heating. Sabalan is a high temperature geothermal region in Ardabil province northwestern of Iran. In this study the conductivity structure of the top crust (depth lower than 2 kilometers) ...
Read More
Geothermal energy is playing a large role as an alternative energy source for both electricity generation and space heating. Sabalan is a high temperature geothermal region in Ardabil province northwestern of Iran. In this study the conductivity structure of the top crust (depth lower than 2 kilometers) is examined using data from collected 212 magnetotelluric (MT) and time-domain electromagnetic (TDEM) soundings across the Sabalan volcano in 1998. The MT data collected between 1-8192 Hz is of useful quality and provides good control on the surface layers in majority of sites. The MT data were corrected for static shift effect using TDEM data. The TDEM data and MT data were jointly inverted to yield 1D and 2D models. In practice for geothermal investigation, measured MT data are never entirely 1D and diagonal elements of the impedance tensor are always non-zero. However, in many cases the Earth response is dominated by an overall 2D structure, permitting data analysis within these lower dimensions. In this study, in order to have the best possible interpretation we used two modes MT data in 2D inversion. Static shift correction and inversion approach accomplished by using WinGlink software. We choose three intersect profiles with 54 MT sites along with, to show resistivity distribution around Sabalan and try to predict the main intrusive magma chamber position as the heat source of the system. MT resistivity images confirmed the findings of previous surveys and the reported geological features in the Sabalan field. To have a better view about resistivity distribution, we combine the 2D models at the intersection points to obtain a 2.5D view about the resistivity in the area. The resulting models reveal the extension of the high conductivity anomalies in the western and southwestern parts of the area that most probably is related to the main heat source of the geothermal system at shallow depths.
M. Montahaei; B. Oskooi; H. Brasse; G. Kapinos
Abstract
Investigations made in current study illustrate that an application of 2-D isotropic inversion algorithm for magnetotelluric data affected by anisotropy could recover macro-anisotropy reasonably well. First anisotropy effects on common MT interpretation steps, dimensionality analysis and 2-D inversion ...
Read More
Investigations made in current study illustrate that an application of 2-D isotropic inversion algorithm for magnetotelluric data affected by anisotropy could recover macro-anisotropy reasonably well. First anisotropy effects on common MT interpretation steps, dimensionality analysis and 2-D inversion modeling, are investigated. Two kinds of azimuthally anisotropic features (anisotropic block and anisotropic layer) which generally form a part of 2-D models are considered. The influence of different anisotropy strikes and resistivity contrasts on dimensionality analysis and on the behavior of induction arrows is studied. These investigations evince that, a strike direction close to the anisotropy strike can be chosen by the dimensionality analysis of the data. Then if the data are rotated to this angle, 2-D inversion would recover the anisotropy sensibly by means of macro-anisotropy. This procedure is tested successfully on a field data set where anisotropy had been previously recognized. The results show that the proposed approach reproduces the anisotropy acceptably via macro-anisotropy in the final inversion model.